Universal Scaling Relations in Strongly Anisotropic Materials 
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We consider the critical temperature in strongly anisotropic antiferromagnetic materials, with 
weak coupling between stacked planes, in order to determine the interplane coupling constant 
from experimentally measured susceptibilities. We present theoretical arguments for a univer- 
sal relation between interplane coupling and susceptibility shown numerically by Yasuda et. al., 
Phys. Rev. Lett. 94, 217201 (2005). We predict a more general scaling function if the system is 
close to a quantum critical point, a similar relation for other susceptibilities than considered in 
Yasuda et. al., and the validity of these relations for more general phase transitions. 
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Many materials display at low temperatures strongly 
spatially anisotropic responses to magnetic or electronic 
probes. This fact has motivated the theoretical study 
of low dimensional quantum systems on their own right. 
Solving one- or two-dimensional quantum systems can 
be useful to understand intermediary regimes of temper- 
ature in which fluctuations are dominated by the sub- 
system of lower dimensionality. Three-dimensionality is 
effectively restored once the temperature is lowered be- 
low the lowest energy scale characterizing the anisotropy. 

For magnetic systems this scale can be the temper- 
ature {//3 AF (k B = h = 1) below which antiferromag- 
netic (AF) long-range order manifests itself, fn this con- 
text, one of the most studied model is perhaps a stacking 
in three dimensions of chains or square lattices on each 
of which a nearest-neighbor quantum spin-S* Heisenberg 
model Hj with AF exchange coupling J > is defined. 
To model a strong spatial anisotropy, one assumes that 
there exists a nearest-neighbor AF exchange coupling 
J' in the directions transverse to the chains or planes 
that is much weaker than J, J 3> J' > 0. The three- 
dimensional quantum Hamiltonian is H 3d . Many efforts 
have been invested for the last 30 years in calculating 
the J'-dependence of 1 //3 AF 0, H H . In this letter 
we shall address the question: are there some universal 
relations that relate J' and some observables of Hj or 

The motivation for this question comes from the work 
by Yasuda et al. in Ref. in which the Neel temperature 
l/(3 AF of H 3d , the n-dimensional static staggered suscep- 
tibility X i n) = X ZZ (Q, u = 0; 8 AF ) of Hj where Q = tt if 
n = 1 or Q = (tt, tt) if n = 2, and J'xs"" 1 = l/C n (J') were 
computed numerically as a function of < J'/ J < 1. In 
quasi-two-dimension, it was observed that 



J'x 



(n=2) 



(1) 



becomes independent of J'/J when J'/J < 0.1 and that 
the constant value that it takes, although not 1/2 as pre- 
dicted by mean-field theory, 0,0 is independent of the 
magnitude of the spin and even takes the same value for 



a classical (S = oo) model. Although the evidence is less 
pronounced the same conclusion was reached in quasi- 
one-dimension. 

We want to construct a tractable model that repro- 
duces qualitatively these findings and we want to under- 
stand how these results can be useful to establish experi- 
mentally the implied universality. We present a theoret- 
ical argument that as J'/J j 0, the function ( n=2 con- 
verges to a constant. We make the following additional 
predictions. 

First, we consider more general AF models in the 
plane, and we consider the case in which, by tuning pa- 
rameters, it is possible to tune the planar model close to a 
quantum phase transition, so that the zero-temperature 
AF order of the two-dimensional model (without inter- 
plane couplings) becomes small. Then, we predict the 
scaling function 



^^(c/W^) 



(2) 



for some scaling function F ± , in the limit J'/J j 0, where 
£( 2 ) is the correlation length of the two-dimensional 
model at temperature \/(3 AF and c is a spin- wave ve- 
locity defined below. Note that in the system considered 
by Yasuda et. al. the planar model is in the renormalized 
classical regime so that cf3 AF /£}- 2 } is exponentially small 
in cj3 AF and converges to zero as J'/J j 0. Therefore, 

^i(c/W£ (2) = o) = i/c ? . _ 

Second, we predict a similar scaling relation that will 
be valid for quantities which are easier to access exper- 

(2) 

imentally. The susceptibility \ s ' defined above is that 
of the two-dimensional model without the interlayer cou- 
plings, and cannot be measured in most real materials. 
We define Xtt.^o = X zz ( 7T , 7r :^,uj = O;0 AF ) to be the 
static susceptibility in the layered system at wave vector 
(tt, tt) in the plane and wave vector perpendicular to 
the plane at temperature l//3 AF . Then, we predict that 



/'v (3) - F (r^B /f^> i 

J A7r,7r,0 — r 2\ c P AF I W,7r,oA 



(3), 
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for some scaling function F 2 , in the limit J'/ J | 0, where 
and are the in-plane spin- wave velocity and cor- 
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relation length of H 3d at temperature l//3 A p near wave 
vector (w, tt, 0), respectively. As it is the instantaneous 
structure factor S n w0 that is most readily measured 0, 
we note that the product J'S n ^ (3af also should obey a 
scaling law of the form © for some scaling function F 3 . 

Third, we predict that similar scaling results hold for 
other layered models. 

In quasi one-dimension, we expect that similar scal- 
ing results will also hold. This does not, however, help 
us understand the results of Yasuda et. al. in quasi one- 
dimension. The scaling functions F 1 , F 2 imply that the 
classical and quantum models will show the same Q- x > 
only if cf3 AF <C £W . However, as the one-dimensional 
Heisenberg model on a chain is not in the renormalized 
classical regime but rather quantum critical, it should 
have some non-zero c0af/£ anc ^ should show a dif- 
ferent C than the classical model. Thus, the one- 
dimensional results remain a puzzle. 

Physical Motivation — Here, we present a physical mo- 
tivation for the results above and a brief microscopic 
derivation of the relevant non-linear sigma model. In the 
next section, we show these scaling results using a renor- 
malization group (RG) for this non-linear sigma model. 
The reason for which we use this model is that we want 
to illustrate the effects of field renormalization and the 
non-linear sigma model RG already has a nontrivial field 
renormalization at leading order in the coupling constant, 
while such a renormalization is not seen until order e 2 
(1/iV) in a 4 — e (large N) expansion. 

Since the interplane interaction is weak, we can treat 
it perturbatively at the microscopic level. Following 
standard steps, in the absence of the interplane inter- 
action, we can first derive the partition function for the 
two-dimensional O(N) quantum non-linear sigma model 
(2e?QNLSM) with field n k (r,T), where k is a discrete 
index labelling individual planes, r is a two-dimensional 
vector describing coordinates in the plane, and r is imag- 
inary time. The relevant action for plane k is S k — 



k — 



(It 



(4a) 



and 



;-(2) 



c 



(2) 
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T Z h h n h . 
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Here, the lattice spacing a plays the role of the micro- 
scopic ultra-violet (UV) cutoff, i.e., A ~ 1/a that of an 
upper cutoff on momenta. The linear size L of the plane 
is the largest length scale of the problem. The deriva- 
tive = (d CT , V) depends on the spin-wave velocity 
c in the plane and is of order J a. The dimensionless 
coupling constant g depends on the microscopic details 



of the intraplane interactions. The dimensionless back- 
ground field h 7 where h = \h\, is the external source for 
a static staggered magnetic field conjugate to the planar 
AF order parameter of the underlying lattice model. It 
breaks the O(N) symmetry of Lagrangian l|4a|l down to 
0(N — 1) and as such acts as an infra-red (IR) regulator. 
The dimensionless coupling Z h is the field renormaliza- 
tion constant associated to n fe . The use of the continuum 
limit within each of the planes labelled by k is justified 
if we are after the physics on length scales much longer 
than a. 

The interplane nearest-neighbor AF coupling J' gives 
the characteristic interplane spin-wave velocity d ~ J'a 
and length scale a' = [J/ J') l l 2 a. The couplings J',g 
get renormalized as discussed below, so the velocity d 
changes at longer length scales. For a very weak nearest- 
neighbor interplain AF coupling, J' <C J, the physics on 
length scales much larger than a but yet not much larger 
than a' is captured by the partition function 



Z-- 



l[V[n k ]S(nl-i: 



(1) 
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The Lagrangian C k encodes the effect of the microscopic 
nearest- neighbor interplane AF interaction J'. To com- 
pute this, we use a Hubbard-Stratonovich transformation 
to replace the microscopic interaction between any two 
spins in the cubic lattice with coordinates (i,j,k) and 
(i,j,k + 1) by an interaction of each spin with a fluc- 
tuating magnetic field H i ^ k . When J' <C J and on 
intermediary length scales a <C A < a', the dominant 
mode for the Hubbard-Stratonovich field is near momen- 
tum (tt, tt) in the plane. The action for the spins in the 
presence of this field is the same as Eq. (|4b|) with field h 
replaced by H k . In this approximation, integrating the 
Hubbard-Stratonovich field gives the short-range inter- 
plane interaction term 



(2) 



c 



(3) 



(5b) 



C 
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J'Z' , N 2 

(n t - n h 
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where Z' renormalizes as Z\ to lowest order in J'/ J, 

Z' = Z 2 h [l + 0(J'/J)}. (6) 

This defines the so-called three-dimensional strongly 
anisotropic O(N) QNLSM (3dSAQNLSM). 

This derivation of the 3g?SAQNLSM considered inter- 
actions between spins positioned directly above and be- 
low each other in neighboring planes. It is possible to 
treat more complicated interplane interactions. For ex- 
ample, going back to the cubic lattice, let there be an 
AF interaction J 1 between the spin at site (i,j,k) with 
that at fc± 1) and another AF interactions J 2 to the 
spins at sites (i ± k ± 1) and (i, j ± 1, k ± 1). Then, 
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if J x < J and J 2 < J, we can derive the 3dSAQNLSM 
with the interplane coupling J' = J 1 — 4 J 2 . 

Renormalization Group — Here, we present an RG 
analysis of the non- linear sigma model (JSJ). The most 
important result in this section is that the identity © 
is preserved under the RG flow up to the length scale at 
which the scale dependent effective anisotropy (|f 3|) is of 
order 1. 

We perform a RG analysis following Polyakov for con- 
venience Q. In each plane labelled by k, we write 





N-l 



n. 



(7) 



a=l 



The field of unit length m k encodes the planar AF or- 
der expected in the limit g/c J, 0, while the N — 1 fields 
e k capture the deviations away from the direction m k 
of symmetry breaking, i.e., the JV-1 ej form an or- 
thonormal basis of vectors orthogonal to m k . The N — 1 
coefficients <f) k make up the vector <p k . To leading order 
in an expansion in powers of g/c of the parametrization 
, the field m k is the slow mode while the N — 1 fields 
<p k represent fast modes with characteristic 2-momenta 
A < \p\ < A. Substituting Eq. J7J) into Eq. © gives the 
Lagrangian 



4 1J 



4 2) 



r>a Tib 



C 
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(8) 



<^<5 Qb )] --Z h h-m fc (l-^) 
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to leading order in an expansion in powers of g/c. 
The N — 1 coefficients are defined by d ll m k — 

T,a=i B ^ e k- The ( N ~ X )( N ~ 2 )/ 2 independent co- 



efficients Af^ 



-(d e%) ■ e%. The RG flows of the di- 



mensionless couplings g, Z h h, and t = 1/{J(3) that follow 
after integration over the fast modes <p k in the limit of no 
interplane interactions were computed by Chakravarty, 
Halperin, and Nelson to leading order in g/c (see Fig.QJ 
||. To this order, c is unchanged. 

To quantify the very weak microscopic interplanar cou- 
pling, we define the anisotropy a as the ratio of the im- 
portance of yCl to C k 3 ^ when the upper cutoff on the 
momenta is A. By assumption, this anisotropy is strong 
at the microscopic level (upper cutoff A), 



a = gJ'Z'/J<€ I. 



(9) 



Next, we consider the renormalization of C k ^ in Eq. @ 
when h = and of the interplane interaction 



r(3) _ 

*^ u — 



J'Z' 
~2aJ 



'l J.2^/2 

1 - 4> k ) m k 



.1/2 



k+lj 



m 



k+l 



N-l 



^fe e fe — 4>k+l e k+lj 



(10) 



FIG. 1: Phase diagram and RG flows for the 3dQNLSM 
(a = 1) and 2dQNLAM (a = 0) after Ref. 0. Conjectured 
phase diagram and RG flow for the anisotropic 3dQNLSM 
(1 > a > 0). The shaded regions have long range Neel order. 
The figure for 1 > a > shows a two-dimensional slice of the 
three-dimensional RG flow, as 5 is changing under this flow. 



after averaging over the fast modes <p k . To this end, we 
shall introduce the renormalized values 



1 

9 
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t 



9 \A 



l + («-<^0 , (lla) 



<Pl6 ab )\ (lib) 



(lie) 



at the scale A as a result of averaging over the fast modes 
<p k . For a <C 1, this average over fast modes is {<j) k <t>\) = 
5 kl S ab ln(A/A)^ coth(#/2i). Furthermore, 



,(3) 



2«2 zi 



m. 



rn 



fe+ly 



from which follows the rcnormalizations 



gZ>_ 

gz r 



z> 
z 7 



7 2 
Zh. 



(12) 



(13) 



so long as a <C 1. 

Let us follow the RG flows encoded by Eqs. i|ll|) 
and Eqs. i|13|) starting from the initial values g > 0, 
t = 1/{J(3 AF ) < oo, and 1 3> a > 0, see Eq. JjJJ, corre- 
sponding to a point on the phase boundary between the 
Neel and paramagnetic phase (see Fig.^l. Aside from the 
thermal de Broglie wavelength of the spin waves c(3 AF , 
the initial values g and t define a second characteristic 
length scale, the correlation length £( 2 ) in the 2dQNLSM, 
in view of 1 ^> a > 0. We shall distinguish two cases. In 
the renormalized classical regime 

quantum critical regime cf3 AF /^ 2 ^ ~ 1. Finally, we de- 
note by £ cross the RG length scale at which a ~ 1 and be- 
yond which the RG flows Eqs. ifTTJl and Eqs. fT3^l should 
be replaced by the flows of the isotropic 3g?QNLSM; 
naive scaling gives Across ~ a', but the RG flows above 
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will change this scaling. Any two of these characteris- 
tic length scales, cpAF^C , and Across, fix the third one 
since the RG flows are constrained to the boundary be- 
tween the Neel and paramagnetic phases by assumption. 
Without loss of generality, we shall consider the case 
A > \/{c(3 AF ). As we lower the upper momentum cut- 
off, the RG scale A -1 will eventually become larger than 
c(3 AF . We shall consider RG scales A -1 » cf3 AF , for 
which the quantum fluctuations are important. 

We begin with the renormalized classical regime of the 
3dSAQNLSM. As is illustrated in Fig. □ the running 
coupling constants g in Eq. (|lf af) and t in Eq. (|llb|l flow 
towards zero and oo, respectively, as A decreases but so 
long A" 1 < cf3 AF . By Eq. (| 1 31) . the effective anisotropy 
decreases; using naive scaling which is valid for g <C f , we 
have a ~ (A/A) 2 a. Beyond the RG length scale A -1 ~ 
cf3 AF we can replace the 2dQNLSM in each plane by a 
classical 2dNLSM with the effective coupling g cl , where 
g cl = g at the scale A -1 ~ cf3 AF . The effective anisotropy 
of the classical 3dSANLSM continues to decrease as 
5 = (A/A)^Ac(3 AF )[g cl Z'/(gZ')]a = (k/ kfg cl Z> J' $ AF 
continues to grow until it reaches the isotropic RG scale 
a *~ 1. Equation (jfJJl can then no longer hold. However, 
since there is only a finite range of scales over which 
a is non-negligible but still less than unity, we deduce 
that, at the scale a ~ 1, Z' ~ Z 2 , up to some constant 
of order unity. Furthermore, g cl ~ 1 at this scale also 
since it lies at some point on the phase boundary be- 
tween the Neel and paramagnetic phases. But g cl ~ 1 
tells us that the corresponding RG scale A is of the order 
of the correlation length of the 2dQNLSM. In turn, this 
allows us to infer that the static staggered susceptibility 
of the 2dQNLSM is given by ~ (A/A) 2 ^/^, up 
to universal corrections of order unity. We now multiply 
Xs by J' estimated from the anisotropy a of the classi- 
cal 3dSANLSM, J' x i 2) ~ a{z\fz')fg d . Using Eq. ®, 
and the fact that g ci ~ I and 5 ~ 1, we arrive at Eq. Q 
in the renormalized classical regime. 

Note that each of these relations, such as a ~ I and 
g cl ~ 1, is defined up to a multiplicative constant that 
depends on the details of how we define the RG. How- 
ever, the dimensionless combination in Eq. Q is uni- 
versal. The reason for the universality is that all the 
microscopic details of the Heisenberg moded are encoded 
into the three independent quantities g, Z', and Z h on 
any length scale much larger than a. Let us perform 
the RG flow to some scale such that a is much less than 
unity. Then, the identity JJjJl relates Z' to Z h , leaving 
only two quantities independent in the classical regime, 
say g cl and a. The requirement of criticality relates g cl to 
a, leaving only one independent quantity, say a. Choos- 
ing the rcnormalization scale to be some given fraction of 
the correlation length in the two-dimensional model fixes 
the last quantity, and thus there are no independent pa- 
rameters left. 



Near a quantum critical point and as is the case for 
the renormalized classical regime, the length scale at 
which a ~ 1 is of the order £' 2 ) . Now, however, there 
is no significant separation of scales between c/3 AF and 
£( 2 ) anymore, i.e., a ~ 1 already at c(3 AF . Correspond- 
ingly, there will be universal corrections to I© in the form 
Z'/Z' = «(1, g{\))Z\jZ\ where the function k of a and 
g is universal with k(0,<?) = 1. The deviations in the 

(2) 

quantum critical regime from the limiting value of J'xs 
in the classical renormalized regime define the universal 
scaling function F 1 of c(3 AF /£}- 2 \ 

Similarly, the correlation at the (it, 7t, 0) point is of 
order A -1 in the plane while it is of order a single in- 
terlayer spacing between the planes. Thus, o ~ X» 
and ^l t0 ~ £( 2 ), and so Eq. © follows. 

We close by noting that all arguments presented here 
for a non- linear sigma model with O(N) symmetry ex- 
tend to non-linear sigma models defined on Ricmannian 
manifolds with a positive curvature tensor. For example, 
we expect similar universal scaling relations for a stack- 
ing of AF Heisenberg models on a triangular lattice. 

Discussion — We have provided a field-theoretic basis 
for understanding the result of Yasuda et. al in the quasi- 
two-dimensional case, generalized it to deal with situa- 
tions near quantum critical points, and given a version 
thereof expressed in terms of experimentally accessible 
quantities. The calculation of the scaling functions F 1 2 3 
within a field theoretic approach will only be approxi- 
mate, and the best estimate for F 1 (0) is given by numer- 
ical calculations. For example, it can be shown that the 
mean-field result F 2 (Q) = 1/2 follows from the large N 
limit of the 0(N) 3dSAQNLSM. It would be very valu- 
able to perform a numerical study of F 2 3 , in both the 
renormalized classical and quantum critical regimes. 
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